Quantum Random Walks of Waves 
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The extremely fascinating behaviors of the quantum walks of particles, which differ much from 
the classical counterparts, have attracted many attentions. Here we investigate another interesting 
part of the quantum walks, that is the quantum walks of waves. Firstly, we show the behaviors of 
the quantum walks of plane wave, which are largely different from the counterparts of either the 
classical or the quantum walks of particles. Two situations - with and without intermediate time 
measurements of the walks are considered. At last, it is shown that the quantum walks of plane 
wave can be used to calculate the evolution of the general wave packets, e.g., Gaussian wave packet. 

PACS numbers: 05.40.Fb, 03.65.Ta, 03.67.-a 



I. INTRODUCTION 

In recent years quantum walks have gained great inter- 
est from physicists, mathematicians, computer scientists, 
and engineers. This is caused by the fascinating prospect 
of the applications of the quantum walks to the algorithm 
in the quantum computer, which is supposed to be born 
in the near future. Recent demonstration of a search al- 
gorithm [l], 0| based on quantum random walks showed 
that it is exponentially faster to solve a problem by the 
quantum walks than by the best classical algorithm. 

The behavior of quantum walks is extremely different 
from the classical random walks. As same as the clas- 
sical case, the quantum walks also have a coin(or many 
coins) and a particle. The particle moves in one direc- 
tion based on the result of the tossed coin. However, in 
quantum case, there can be no measurement of the coin 
states during the process of the walks, moreover, both the 
flip of the coin and the conditional motion of the particle 
are unitary transformations, so the process is reversible. 
The variance of the probability of the particle's position, 
which grows linearly with time in the classical walk, by 
contrast, grows quadratically with time in the quantum 
walk, for the reason of the interference between the pos- 
sible paths of the particle in the quantum walks when 
there is no measurement at intermediate time steps. 

The quantum walks of particles have been studied very 
well in recent years. Since the quantum random walks 
has been proposed by Aharonov et al. [3|, many as- 
pects of the quantum walks have been much studied. 
Nayak and Vishwanath have got the exact solutions of 
the Hadamard walk of particles in one dimension in 
Ref. Quantum walks include one dimension 0-[ll| 
and high dimensions in graphs (L2hT3 |. discrete and con- 
tinuous [lij], periodic and quasi-periodic or random 
etc. T. Brun et al. have studied much on the multi- 
coin and decoherence walks 0-0, Qj} , and the absorbing 
problems have been well discussed in pjj [3| ■ 
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However, the quantum walks of waves have not been 
much studied yet. Using the long-wave approximation, 
Aharonov et al. have studied the random walks of a 
Gaussian wave packet in Ref. @ , in which the walk is set 
to be measured after each step. We will review this work 
in the following sections. 

Here, we proceed as follows: first we give a brief re- 
view of the quantum walks of particles in section II, in 
which we also do some improvement of the walks. On the 
one hand, this section would acts as a contrast with the 
quantum walks of waves described in the following sec- 
tions; on the other hand, it would be more convenient to 
describe the walks of waves by introducing the walks of 
particles first. In section III we would turn to the study 
of quantum walks of waves, which consists of the brief 
introduction to the background of the quantum walks of 
waves proposed in Ref. [3| and the study of quantum ran- 
dom walks of plane wave, including both the cases with 
measurements and without any measurement in the in- 
termediate time steps. At last, we will use the results of 
walks of plane wave to study the walks of general wave 
packets in section IV, in which we will give some numer- 
ical simulation of the result and compare some of them 
with Ref. 0. 



II. THE QUANTUM WALKS OF PARTICLES 

In the discrete quantum walks, each step consists of 
two parts: the flip of coin and the condition motion of 
the particle. The value of the coin controls the motion of 
the particle, with the state \R) to go the right direction 
by a distance of I, and the state \L) means to go left. 
Generally, quantum walks are directed by the unitary 
evolution operator. 

For an initial state: 
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the unitary evolution operator can be described as 

E=(P R ®S + P L ®&){C®i). (2) 
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Here the general unitary coin operator 
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The space operator is 
with 

S = e~ lpl/h , S 4 = e ipl/n 
S\x) = \x + l), &\x) = \x-l), 



(3) 
(4) 

(5) 



P R = \R)(R\, P L = \L)(L\. 

After t steps, without any measurements during the 
process of the walks, the state of the particle becomes 

|*(t)> = #|*(0)> = J2[a R (x,t)\R) + a L (x,t)\L)]^ \x). 



The general case with three variable parameters ry, 0, 
and 9 has been studied in the Ref. One special case, 
with </) = rj = 6 = 0, 1 = 1, which is the well known 

Hadamard walk H — ( J ^ \ , has been much stud- 
ied [1, [j| . By either the way of Fourier analysis or the 
way of combinatorial analysis one can get the following 
analytical result: (in Ref. Q) 



|0>® \R) a R (x,t)\R) + a L (x,t)\L), (6) 



|0)<8|L) 



H 



t steps 
H 



>b R {x,t)\R)+b L (x,t)\L). (7) 



However, if we reverse the sequence of the space opera- 
tor and coin operator, that is, we can carry the condition 
motion first and then use the coin flip at each step, an- 
other very different result will be obtained. 

This process may be useful in some cases, for example, 
the measurement of a component of the spin of a spin- 
1/2 particle. One particle moves conditionally first in the 
electromagnetic field by the different force caused by the 
different spin states in the electromagnetic field, and then 
make the spin transformation which is correspondent to 
the coin flip. For the general case of <fi, 77, 9, one can get 
the result in the same way. 



III. THE QUANTUM WALKS OF COSINE 
WAVE 



We start from the introduction of the work Aharonov 
et al. have done in Ref. [3j, which is mathematically 
equivalent to the process introduced in Ref. . There is 



a wave packet with the coin state l^m) and initial space 
state f(x,0). The walk of the wave packet at each step 
consists of two parts: letting the operators U and R{9) 
act on its state sequentially, and measuring both coin 
and space states and then re-initializing the coin state 
into |* in ). 



|* in ) = a R \R)+a L \L), 



U = \R)(R\ ® e~ lpl/h + \L)(L\ <g> e tpl/h , 



R(0) 



cost 
sinf 



— sin 9 
cos 9 



(8) 



(9) 



(10) 



However, to get a good result of the walks, one must 
use the long-wave approximation, that is, the width of 
the wave is much larger than the step length /. By this 
approximation, the wave packet will just be translated in 
space with no variance of shape each step. What's more, 
the moving distance in one step can be much larger than 
the step length I with a low probability if we choose some 
special values of 9. The results of Gaussian wave without 
long-wave approximation can be simulated in numerical 
method. 

Now we are interested in two things. First, for a gen- 
eral wave packet, what the results will be without any ap- 
proximation. Second, what the wave packet will be after 
t steps without measurement in the process of walking. 
However, the direct calculation is very complicated and 
difficult, so we will study the long- wave limit — the plane 
wave first, and then expand the results to the general 
wave packets. (In contrast, the short-wave limit, which 
is actually particles, has been reviewed in section II.) 

Before we study the quantum walks of plane waves, it 
is necessary to note the fact that two plane waves with 
the same frequency can be combined into one plane wave 
with no change in frequency, i.e., 



with 



C = y/A 2 + B 2 + 2AB cos k(a-b), 



tan kc = tan 



A sin ka + B sin kb 
Acoska + B cos kb 



: tan 



arctan(a//3) 
kc = { 7r + arctan(a//3) 
arctan(a//3) — n 



P > 

< 0,a > 
/3 < 0,a < 



(11) 



(12) 



(13) 



(14) 
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A. With measurement after each step 



It is assumed that the initial state of the plane wave 
is |^in) <8> \e lk ( x ~ l °)), then the operate U and operator 
R(9) are carried on it sequentially, the coin state will be 
measured and re-initialized into l^in)- After repeating 
this work several times the state will become l^in) ® 
| e ifcO-2 )y ]\j ow we consider the next step of the wave: 

I*') = U\^in) <8 \e lk{x - V) ) 

= a R e lk{x - l '- l) \R) +a L e tk{x - l ' +l ^\L), 



msm = ( cos , ~ sm ; 

w 1 ' \ smO cosO 



-. -ikix — V — l) 



a L e 



ik(x — l' 



/pie 



ik{x — l ! 

ik{x-l'+l 2 ) I ' 



(15) 



with 

Pr = \ a R cos6*| 2 + \a L sin6>| 2 — 2a R a L sin#cos#cos2fcZ, 

(16) 

Pl = |a H sin6>| 2 + \a L cos#| 2 + 2a R a L sin 9 cos 9 cos 2kl, 

(17) 



— a R costl — a L sin(/ 

tanfcii = : — — tanfet, (18) 

a R cos 9 — a, sin 9 



-a R sin 9 + a L cos 9 

tanK/2 = : — - — tanKf. (19) 

a R sm 9 + a L cos 9 

We can get an amazing result: from the equations above 
we can see that the walking probability to the left or the 
right and the walking length are independent of the steps 
t. That is to say, in the whole process the frequency of 
the wave does not change, and in each step the plane 
wave moves to right or left with a fixed probability, and 
the moving distance is fixed too. This is similar to the 
result achieved in Ref. [3|, but there is no approximation 
here. 

After t steps, the average moving distance of the phase 
of the plane wave is: (x) = t(p R l\ + p^h) , the variance 
cr 2 = tp R p L (li — l 2 ) 2 - By the way, the value of x + 2nir/k 
has the same meaning with the value of x for plane wave. 




FIG. 1: The values of L\ and L2 act as the function of 9. 
With k = 1, I = 0.01, a R = at, = 1/V2. 
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FIG. 2: The values of Li and Li act as the function of fre- 
quency k. With I — 9 — 1, a R — a L = l/\/2. 



and L2 acting as the function of 9. We can see, both 
Li and Li are almost zero for some value of 9, while in 
some other areas one of them would be much larger than 
the other. We choose 9 = 5.55 to plot the relationship 
between L\ L 2 and k in Fig. 2, where L 2 is much larger 
than L\. 



Now, let's discuss the result. We care about the way 
the values of l\ and I2 varying with the variance of k 
and 9, for two reasons. First, one wave packet can be 
expanded in the basis of plane waves with different fre- 
quencies rang from —00 to +00. Second, in some situ- 
ation in physics, one may consider the propagation of a 
wave packet, and in the route of the propagation there are 
some disturbance which is equal to the quantum walks 
introduced above. What's more, the disturbance is vari- 
able or even random, which equals to the variance or 
randomness of the value of 9. So in Fig. 1 we plot L\ 



B. No measurement at intermediate time steps 

Now we improve the walks in Ref. Q to the non- 
measurement case, which means we will not measure the 
state of the plane wave in the intermediate process, and 
let it evolute with the interference between the possible 
paths. This result would be more complicated. 

We suppose the initial state is 

|*(0)) = (a R \R) + a L \L)) ® e ^ x ~ l °\ (20) 
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the space operator is 

U = S ® P R + 5 f ® P L 

= |i?)(i?|®e-^ /R + |£)(L|®e^ /ri 







iPl/h 



(21) 



and the evolution operator is 
E' = R(9)S = 



cos f — sm 
sin 9 cos 




,-iPl/H 





jpi/h 



-iPl/h CQS Q _ e iPi/h i 



I e -iPl/h sin $ ^Pi/ft 
The state after i steps will be 



cos t 



!*(*)> = # |*(o)) 



(22) 



(23) 



To solve this equation we can get the eigenstates and 
eigenvalues of matrix M first: 



M 



= Ai 



ill 



= A 2 



1^(0)) can be expanded in the two eigenstates: 



p ik(x—l ) 

so the final state is 



m)) = E' t \m) = (l R L ( § ) )=v> R (t)\R)+v>,(t)\L) 



= A l ( £ ) 771 + A 2 



Through much calculation we can get the analytical 
solutions: 

where the value of Pi, P2, L\, and L2 are determined 
by 9, k and the steps t. The details of the mathematical 
expressions and derivations of them are in the Appendix. 

The result shows that in the end of the walking the 
plane wave is still a plane wave, with a change in phase 
and no variance in frequency. We will use these results 
to study the quantum walks of general waves. 



IV. THE QUANTUM WALK OF GENERAL 
WAVE PACKETS 



As we know, every wave packet can be expanded in the 
form of plane waves (or equivalently cosine wave): 



/(z,0) 



+ OO 



f(k)e lkx dk 



(25) 



where f(x, 0) is the initial space state of the wave packet. 
As long as we know the evolution of the plane wave e lkx , 
we can use the integral (25) to get the evolution result of 
wave f(x), and the f(k) can be treated as the weight of 
the plane wave e lkx . 



1 



+ OO 



/( fc ) = ^/ f(x,0)e~ ikx dx 



(26) 



A. Do measuring after each step 

After one step, the plane wave e lkx develops into: 

y/p^R) ® e lk(x+h) + y/p^L) ® e lk{x+h) . (27) 

After t steps, with measuring the coin state each step, 
the wave packet would develop into 



with the probability p 



e ik(x+nh+(t-n)l 2 ) ) n g [0, t] 
t 



(28) 



n 



p\p L ™ . So one can study 



the quantum walks of any wave packet with the interme- 
diate measurements, using Eq. (25). 




FIG. 3: The probability amplitude distribution for the Gaus- 
sian wave, with the same condition in the Ref. [j§], we choose 
the parameters 6 = — arctan(0.9), I = 0.01 to make the dis- 
tance the wave moves each step much larger than I. This 
figure gives the shapes of the packet after 1, 3, 5, 10, and 20 
steps, respectively. 

For a clear insight, we can give a typical example 
- the walks of the Gaussian wave packet f(x,0) = 
cxp (— x 2 /2) /7T 1 / 4 . We consider the case that the re- 
sult of the coin state measured each step would be \L). 
That is to say, the wave moves left each step. If the wave 
moves left each step, then the result after t steps will be 

r+00 

f(x, 0)^C f(k) (v?!)' \L) ® e lk(x+th Uk , (29) 
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FIG. 4: The probability amplitude distribution for the Gaus- 
sian wave after 35 steps. 

where C is the normalization coefficient, to make sure 
the whole probability of the wave packet be unit; p L and 
?2 are given by equations (17) and (19); and by Eq. (26), 
we have 

f(k) = l/V2^cxp(-k 2 /2). (30) 

Fig. 3 and Fig. 4 indicate that the wave will be split into 
several small waves after many steps of quantum walks. 
If 9 is random and changing at each step, the results will 
be more interesting, here we will not proceed any more. 

B. No measurement at intermediate time steps 

Assume that the initial state is 

{a R \R)+a L \L))®f{x,0), 
after t steps it will evolute into 

j f(k) (v^i?) ® e ik ^ x+L ^ + y/ft\L) <g> e ifc ( x+ia >) dfc, 

— oo 

(31) 

where the parameters Pi, Li, P2, L 2 are functions of k 
and 9, which are given in the Appendix. 

We give the results by simulating the Gaussian wave: 
as same as the former case, we choose the parame- 
ters 9 = — arctan0.9, I = 0.01, a R = a L = l/V%, 
f(x,0) = exp(— x 2 /2)/7r 1 / 4 . In figure 5 we note the fact 
that after one step the result is the same to the case of 
doing measurement at each step. We also find that when 
t is not large enough the wave packet will not change 
much, and then we plotted the evolution of the Gaussian 
wave packet after 6000 steps in Fig. 6 to 9, the shapes of 
waves have much difference with initial shape, and would 
be split to some small waves. These figures indicate that 
when the wave has only one peak, the quantum walks 




x 



FIG. 5: The probability amplitude distribution for the Gaus- 
sian wave after one step. The result is the same to the former 
case. One wave moves left with much larger distance than 
and the other moves right with the distance similar to I and 
so almost coincide with the initial Gaussian wave packet. 




FIG. 6: The probability amplitude distribution for the Gaus- 
sian wave after 6000 to 6003 steps. The solid lines represent 
the initial Gauss wave; the dotted-dashed lines represent the 
waves with \R) coin state; the dashed lines represent the waves 
with \L) coin state. 

will make little change on its shape at each step, but if a 
wave has two or more peaks the quantum walk will make 
large difference to its shape at each step for the reason 
that these peaks may interference each other during the 
walks. 



V. CONCLUSIONS 

We have introduced a new idea — the quantum walks 
of waves, and proposed a method to study and solve this 
problem. First, to simplify the calculation we studied the 
walks of plane wave, in which we considered two situa- 
tions: with measurement in each step and without mea- 



suremcnt in the intermediate process. Hence, we studied 
the quantum walks of general waves using the results of 
plane wave, which also includes the two situation men- 
tioned above. Moreover, we have simulated some special 
wave packets. In the whole process, we find the behav- 
ior of waves quantum walking is much different from the 
counterparts of particles and also classical situation. By 
the way, the short-wave limit of wave is just delta func- 
tion, which is equivalently particle. 

The quantum walks of waves can be used in many ar- 
eas of physics, such as the propagation of light field in 
fibers, the transmission of signal with a width, and so 
on. In the propagation and transmission, there would 
be much disturbance such as reflection, refraction, and 
scattering, which are equal to the quantum walks with 
variable values of 9. That is to say, the value of 9 would 
be random in the process of quantum walking. 
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Ai, 2 = L -(e iPl ' h + e- iPl ' h )cosO 

T l^( e iPi/fi - e -iPi/ft)2 cos 2 _ 4 sin 2 fl^ 



with 



a = arctan ■ 



cos kl cos 9 



(32) 



\J sin 2 kl cos 2 9 + sin 2 9 
After t steps, the wave develops into 

<P R (t) = ^fp[e lk ^- lo+L ^ and ip L (t) = ^/p^e ik ^- lo+L ^ 
If t is even P 1 = A 2 + B 2 , P 2 = C 2 + D 2 , 

tan(fcLi - tw/2) = B/A, tan(fci 2 - tw/2) = D/C. 
If i is odd P 1= E 2 + F 2 , P 2 = G 2 + H 2 , 

tan(fc£i - i7r/2) = F/E, t&n(kL 2 - tn/2) = H/G. 



with 



A = a R costa + a L cos kl sin 6* sinta/Q, 



B = {a R cos 9 + a L sin 9) sinklsmta/Q, 



APPENDIX: THE RESULT OF PLANE WAVE 
WITHOUT MEASUREMENT 



C = a L cos to — a R cos kl sin 9 sin ta/Q, 



To get the eigenvalues and eigenstates of M, we notice 
that, for the plane wave e lkx , the operator e ± lPl / h is 
equivalent to e ±lkl . 



M = 



-iPl/h 



COS! 



. e iPl/h sine 



M , , 

h 

We can get: 



Ai 



M 



A 2 



D = (a R sm9 — a L cos 9) sinfc/sinto/Q, 



E = — {a R cos 9 + a L sin 9) sin kl cos ta/Q, 



F = —a R smta + a L cos kl sin 9 cos ta/Q, 



G = (—a R sin9 + a L cos9) sinkl costa/Q, 



a = e lpl/h sin9, 



H = —a L sinto — a H cosfc/sin#costo/Q. 



_ {e iPl/h_ e -iPl/h )cos0 



± -jy{^ Pllh - e- lpl / h ) 2 cos 2 9 - 4sin 2 9, 



Here 



Q = Vsin 2 kl cos 2 9 + sin 2 9. 
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